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Abstract 

In this paper we discuss a novel procedure for constructing clusters of bound particles 
in the case of a quantum integrable derivative 5-function Bose gas in one dimension. It 
is shown that clusters of bound particles can be constructed for this Bose gas for some 
special values of the coupling constant, by taking the quasi- momenta associated with 
the corresponding Bethe state to be equidistant points on a single circle in the complex 
momentum plane. We also establish a connection between these special values of the 
coupling constant and some fractions belonging to the Farey sequences in number theory. 
This connection leads to a classification of the clusters of bound particles associated with 
the derivative (5-function Bose gas and allows us to study various properties of these 
clusters like their size and their stability under the variation of the coupling constant. 
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1 Introduction 



Exact solutions of one- dimensional (ID) quantum integrable many-body systems with 
short range interactions have emerged as an active area of research PH2U]. due to their 
effectiveness in describing recent experiments using strongly interacting ultracold atomic 
gases [2T1 - [28] . The nearly ID motion of ultracold bosons is achieved in such experiments 
by confining the bosons within waveguides that tightly trap their motion in the two 
transverse directions, and allow them to move only in the third direction. Many results 
of these experiments can be understood within the framework of the well known Lieb- 
Liniger model or the 5-function Bose gas; this is a ID quantum integrable system with 
the Hamiltonian for N particles given by 

ff^ = -n2^^ + 2nV ^5(x/-x™), (1.1) 

j=l l<m 

where fi is the coupling constant. Exact eigenfunctions of the Hamiltonian (11. ip are 
constructed by using the methods of coordinate as well as algebraic Bethe ansatz [THn]. 
Moreover, the equilibrium properties of this system have been studied by employing 
the thermodynamic Bethe ansatz and various correlation functions have been computed 
through different approaches [^[8]-[20]. 

In this context it may be recalled that, for a class of exactly solvable dynamical 
systems, the method of coordinate Bethe ansatz directly yields the eigenfunctions in the 
coordinate representation. One can study the asymptotic form of these eigenfunctions 
in the limit of infinite length of the system (i.e., when all Xj's are allowed to take 
value in the range — oo < Xi < oo). If the probability density associated with an 
eigenfunction decays sufficiently fast when any of the particle coordinates tends towards 
infinity (keeping the centre of mass coordinate fixed, for a translationally invariant 
system), a bound state is formed. The stability of a bound state, in the presence of 
small external perturbations, can be determined by calculating its binding energy. It 
is well known, for the case of the (5-function Bose gas (11.11) with N > 2, that bound 
states with positive binding energies exist for all negative values of the coupling constant 
/i [2l[3l|5H9]. The quasi-momenta associated with such a bound state are represented 
by equidistant points lying on a straight line or 'string' parallel to the imaginary axis 
in the complex momentum plane. Moreover, for the case of the (5-function Bose gas 
with negative values of the coupling constant, one can construct Bethe eigenfunctions 
corresponding to more complex structures like clusters of bound particles and show 
that those clusters of bound particles are stable under scattering. The quasi-momenta 
corresponding to such clusters of bound particles are represented through discrete points 
lying on several 'strings', all of which are parallel to the imaginary axis in the complex 
momentum plane [3t[9l [T2HT5] . 

Similar to the case of the (5-function Bose gas mentioned above, there exists another 
exactly solvable and quantum integrable bosonic system with a Hamiltonian given by 

j=l J Km 
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where is a real nonzero coupling constant [291433] . The Hamiltonian (11 ■2p of this deriva- 
tive 5-function Bose gas can be obtained by projecting that of an integrable derivative 
nonlinear Schrodinger (DNLS) quantum field model on the A^-particle subspace. Clas- 
sical and quantum versions of such DNLS field models have found applications in differ- 
ent areas of physics like circularly polarized nonlinear Alfven waves in plasma, quantum 
properties of optical solitons in fibers, and in some chiral Tomonaga-Luttinger liquids ob- 
tained from the Chern-Simons model defined in two dimensions [331flT] . The scattering 
and bound states of the derivative 5-function Bose gas fll.2p have been studied exten- 
sively by using the methods of coordinate as well as algebraic Bethe ansatz [2nH5Slll2Hll] • 
It turns out that the quasi-momenta associated with a bound state of this model can 
be represented by equidistant points on a circle or circular 'string' with arbitrary radius 
and having its centre at the origin of the complex momentum plane. For the cases N = 2 
and = 3, this type of bound states can be constructed for any value of rj within its full 
range: <| ?7 |< oo. However, for any given value of > 4, the derivative (5-function 
Bose gas allows bound states in only certain non-overlapping ranges of the coupling con- 
stant rj (the union of these ranges yields a proper subset of the full range of rj), and such 
non-overlapping ranges of f] can be determined by using the Farey sequences in number 
theory |32Hllj. Furthermore, for any given value of A^ > 3, bound states with positive 
as well as negative binding energies can be constructed. From the above discussions it 
is evident that the bound states of this derivative 5-function Bose gas exhibit a much 
richer structure in comparison to the case of the 5-function Bose gas. 

The aim of the present work is to explore how clusters of bound particles can be 
constructed in the simplest possible way for the derivative (5-function Bose gas (II. 2p . In 
analogy with the case of the 5-function Bose gas, one may think that clusters of bound 
particles can only be constructed for the above mentioned case by properly assigning 
the corresponding quasi-momenta on several concentric circles or circular 'strings' in 
the complex momentum plane. However in this article we shall show that for the 
derivative (5-function Bose gas with some special values of the coupling constant r], 
clusters of bound particles can be constructed in a much simpler way by assigning the 
corresponding quasi-momenta as equidistant points on a single circle having its centre 
at the origin of the complex momentum plane. The arrangement of this article is as 
follows. In Sec. 2, we first briefly review the construction of Bethe eigenstates for the 
case of the derivative 5-function Bose gas. Then we consider a sufficient condition for 
which a Bethe eigenstate would represent clusters of bound particles. This sufficient 
condition for obtaining clusters of bound particles has not attracted much attention in 
the literature, probably because it does not yield any solution at all for the case of the S- 
function Bose gas. However, in Sec. 3, we show that this sufficient condition yields many 
nontrivial solutions for the case of the derivative (5-function Bose gas. Subsequently, by 
using some properties of the Farey sequence, we classify all possible solutions of this 
sufficient condition and obtain different types of clusters of bound particles. In Sec. 4 
we discuss various properties of such clusters of bound particles, such as the sizes of the 
clusters, their stability under the variation of the coupling constant, and their binding 
energy. We end with some concluding remarks in Sec. 5. 
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2 Conditions for forming clusters of bound particles 



In the coordinate representation, the eigenvalue equation for the Hamiltonian ( II ■2p of 
the derivative (5-function Bose gas may be written as 



where ti^{xi,X2, ■ ■ ■ ,xn) is a completely symmetric A^-particle wave function. Since T-Ljy 
commutes with the total momentum operator given by 



tn{xi,X2, ■ ■ ■ ,xn) can be chosen as a simultaneous eigenfunction of these two com- 
muting operators. It may be noted that the Hamiltonian (11.21) and momentum (12.21) 
operators enjoy the scaling property Ti^ — )■ X^Hn and Vn — ^ APat, when all the coor- 
dinates are transformed i/X. Hence, from any given eigenfunction of and 
Vn, one can generate a one-parameter family of eigenfunctions by scaling all the Xj. It 
may also be observed that Hn remains invariant while Vn changes sign if we change 
the sign of t] and transform all the — )■ —Xi at the same time; such a transformation 
may be called as 'parity transformation'. Due to the invariance oH-Ln under this parity 
transformation, it is sufficient to study the eigenvalue problem (12.11) for one particular 
sign of 7], say, > 0. The eigenfunctions for < case can then be constructed from 
those for r/ > case by simply changing Xi — —Xi\ this leaves all energy eigenvalues 
invariant but reverses the sign of the corresponding momentum eigenvalues. 

For the purpose of solving the eigenvalue problem (12.11) through the coordinate Bethe 
ansatz, it is convenient to divide the coordinate space i?^ = {xi^X2, ■ ■ ■ xn} into various 
A^- dimensional sectors defined through inequalities like < x^(2) < ••■ < x^^^n), 

where {cj(l),a;(2), ■■■ ,uj{N)} represents a permutation of the integers {1,2, ■■■ ,N}. 
Since the interaction part of the Hamiltonian (II. 2p vanishes within each such sector, 
the resulting eigenfunction can be expressed as a superposition of free particle wave 
functions. The coefficients associated with these free particle wave functions can be 
computed by using the interaction part of the Hamiltonian (11.21) . which is nontrivial only 
at the boundary of two adjacent sectors. It is known that all such coefficients, which 
appear in the Bethe ansatz solution of a A^-particle system having only local interactions 
(like the (5-function or derivative (5-function type interactions), can be obtained by simply 
solving the corresponding two-particle problem |25]. Thus, by using the solutions of 
the related two-particle problem, it is possible to construct completely symmetric N- 
particle eigenfunctions for the Hamiltonian (11.21) . In the region xi < X2 < ■ ■ ■ < xat, 
such eigenfunctions can be written in the form [291I3T] 



T-Ln tn{xi,X2, ■ ■ ■ , Xn) 



E tn{xi,X2, ■■■ ,Xn) 





(2.2) 



Tn{xi,X2, ■■■ ,Xn) = n 



uj \l<m 



A(^km, ki) 



) 



Pw(l),w(2),- ,uo(N){Xi, X2, ■ ■ ■ 



Xn) , (2.3) 



where 




(2.4) 
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knS are all distinct quasi-momenta, 00 represents an element of the permutation group 
for the integers {1,2, ....A^} and implies summing over all such permutations. The 
coefficient A{ki, k^) in Eq. fl2.3l) is obtained by solving the two-particle problem related 
to the derivative 5-function Bose gas and this coefficient is given by 

The eigenvalues of the momentum (12. 2p and Hamiltonian (II. 2p operators, corresponding 
to the eigenfunctions r7v(a;i,X2, ■ ■ ■ ,xn) of the form (12. 3p . are easily obtained as 

N 

Vn TNixi,X2,- ■ ■ ,xn) = h(^^kj^ tn{xi,X2,- ■ ■ ,xn) , (2.6a) 

i=i 

N 

Hn tn{xi,X2,- ■ ■ ,xn) = h'^(^^ k"^^ tn {xi,X2,- ■■ ,xn) ■ (2.66) 

i=i 



It should be noted that, Bethe states of the form (12. 3p represent scattering as well as 
bound states for the Hamiltonian (II. 2p of the derivative 5-function Bose gas. However, 
for the case of scattering states all the fcj's are real numbers, while for the case of bound 
states the kj^s are allowed to take complex values in general. As mentioned above, 
for a translationally invariant system, a wave function represents a localized bound 
state if the corresponding probability density decays sufficiently fast when any of the 
relative coordinates measuring the distance between a pair of particles tends towards 
infinity. To obtain the condition for which the Bethe state (12. 3p would represent such 
a localized bound state, let us first consider the following wave function in the region 
Xi < X2 < ■ ■ ■ < Xn '■ 

N 

Pi,2,-,n{ xi,X2,- ■ ■ ,xn ) = exp{i^kjXj), (2.7) 

i=i 

where kj^s in general are complex valued wave numbers. As before, the momentum 
eigenvalue corresponding to this wave function is given by ^Xl^Li^j- Since this must 
be a real quantity, one obtains the condition 

N 

E^i = 0, (2.8) 

where qj denotes the imaginary part of kj. By using the condition (12. 8p . the probability 
density corresponding to the wave function pi,2,- -,Ar( Xi,X2,-- - ,XN)m (12. 7p can be 
expressed as 

A^-l r 

\Pi,2,-,n{xi,X2,--- ,XN)f = exp I 2^ (^^gjj ?/r I , (2.9) 

r=l j=l 
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where the j/r-'s are the — 1 relative coordinates: yr = Xr+i — Xr ■ It is evident that the 
probabihty density in fl2.9p decays exponentially in the limit i/^ — t- oo for one or more 
values of r, provided that all the following conditions are satisfied: 



Af-l 



gi<0, gi + g2<0, , < . (2.10) 

i=i 

It should be observed that the wave function (12.71) is obtained by taking u as the identity 
permutation in (12. 4p . However, the Bethe state (12. 3 p also contains terms like (12. 4p with 
u representing all possible nontrivial permutations. The conditions which ensure the 
decay of such a term, associated with any nontrivial permutation a;, are evidently given 

by 

N-l 

qu{i) < , ga;(i) + quj{2) < , , ^ g<^(j) < . (2.11) 

i=i 

It is easy to check that above conditions, in general, contradict the conditions given in 
Eq. (I2.10p . To bypass this problem and ensure an overall decaying wave function (12. 3p . 
it is sufficient to assume that the coefficients of all terms Puj{i),u{2),--- ,ijj{N){.xi,X2, ■ ■ ■ , xn) 
with nontrivial permutations take the zero value. This leads to a set of relations given 

by 

A{K,kr+i) = 0, for reVtN, (2.12) 

where VLn = {1, 2, ■■■ , A^ — 1}. Thus the simultaneous validity of the conditions (12. Sp . 
(12.1 op and (12.120 ensures that the Bethe state tn{xi^X2, ■ ■ ■ ,xn) (12. Sp represents a 
bound state. 

Let us now discuss how the conditions (12. 8p . (I2.10p and (12.120 can be simplified for 
the case of the derivative 5-function Bose gas. Using the conditions (12.80 and (12.120 
along with Eq. (12.50 . one can easily derive an expression for all the quasi-momenta as 

= X e-^(^+i-2")^ , (2.13) 

where x is a real, non-zero parameter, and (p is related to the coupling constant t] as 

(j) = tan~^(?7) =^ ?7 = tan0. (2.14) 

To obtain an unique value of from the above equation, it may be restricted to the 
fundamental region — | < 0(7^ 0) < |. Furthermore, since we have seen that Hn 
(11.20 remains invariant under the 'parity transformation', it is enough to study the 
corresponding eigenvalue problem only within the range < < |. Next, let us 
consider the remaining conditions (I2.10p for the existence of a localized bound state. 
Since summation over the imaginary parts of fc^'s in (I2.13P yields 

E^. = -X^ sin[(iV-O0], (2.15) 

one can rewrite the conditions ( 12.10p in the form 

X ^^^^ sin[(A^-O0] > 0, for / G fi^ . (2.16) 
sin w 
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where Qn denotes the set of integers {1,2, ■■■,A^ — 1}. Consequently, for any given 
values of and A^, a bound state will exist when all the inequalities in Eq. f l2.16p are 
simultaneously satisfied for some real non-zero value of x- 

Let us make a comment at this point. Within the region xi < X2 < ■ ■ ■ < X]\f, the 
completely symmetric Bethe eigenfunctions associated with the Hamiltonian ( II. ip of the 
5-function Bose gas can also be expressed through Eq. (12. 3p . where A{ki,km) is given 



Hence, the equations (12.81) . (I2.10p and (I2.12p together give a sufficient condition for the 
existence of bound states for the case of the 5-function Bose gas also. By using the 
above mentioned form of A{ki, k^), it is easy to show that Eq. (I2.12p completely fixes 
the corresponding quasi-momenta as equidistant points on a straight line parallel to the 
imaginary axis in the complex momentum plane. Moreover, the condition (12. 8 p ensures 
that such equidistant points would be symmetric under refiection with respect to the real 
axis in the complex momentum plane. It is easy to check that the remaining conditions 
(I2.10p for bound state formation are trivially satisfied by those quasi-momenta for any 
negative value of the coupling constant fi. Thus, for any N > 2, bound states are formed 
for the case of the 5-function Bose gas for all negative values of the coupling constant. 
However in the previous paragraph we have seen that, for the case of the derivative 
5-function Bose gas, the conditions (12. 8 p and (I2.12p lead to Eq. (I2.13P implying that 
the quasi-momenta associated with a bound state are represented through equidistant 
points on a circle or circular 'string' with an arbitrary radius and having its centre at 
the origin of the complex momentum plane. Consequently, the remaining conditions 
(I2.10p for bound state formation yield Eq. (I2.16p . which is quite nontrivial in nature. 
Indeed, by solving this equation, we have found earlier that the derivative (5-function 
Bose gas allows bound states in only certain non-overlapping ranges of the coupling 
constant r) Furthermore, such non-overlapping ranges of t] crucially depend on 

the value of and they can be determined by using the Farey sequences in number 
theory. 

We would now like to find the conditions for constructing clusters of bound particles 
in the case of the derivative 5-function Bose gas. To this end, we shall first discuss the 
concept of a 'clustered state' for any translationally invariant system, and then give a 
prescription for finding the Bethe states representing clusters of bound particles. Let us 
consider a system of particles which are divided into some groups or clusters — with 
at least one group containing more than one particle. It is assumed that particles within 
the same group behave like the constituents of a bound state, but particles corresponding 
to different groups behave like the constituents of a scattering state. More precisely, a 
wave function corresponding to such an iV-particle system satisfies the following two 
conditions. If the relative distance between any two particles belonging to the same 
group goes to infinity, the probability density corresponding to the A^-particle wave 
function decays in the same way as a bound state. On the other hand, if the relative 
distance between any two particles belonging to different groups tends towards infinity, 
the probability density remains finite similar to a scattering state. If any wave function 
corresponding to a A^-particle system satisfies these two conditions, we define it as a 



by [I1E1I29] 
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clustered state. It is interesting to observe that, for the case of the 5-function Bose gas, 
the Bethe states corresponding to clusters of bound particles can always be expressed as 
linear superpositions of several clustered and bound states, with the restriction that at 
least one clustered state must be present in such a superposition piTlj . This observation 
may be used to define the Bethe states corresponding to clusters of bound particles for 
the case of any translationally invariant system. More precisely, we assume that a Bethe 
state for any translationally invariant system would represent clusters of bound particles 
if it can be expressed as some linear superposition of clustered and bound states, with 
the restriction that at least one clustered state must be present in this superposition. 

Next, let us discuss how the conditions for constructing a clustered state can be 
implemented for the case of the plane wave function (12. 7p . Since any eigenvalue of the 
momentum operator Vn must be a real quantity, Eq. (12. 8p is also obeyed for this case. 
To proceed further, let us choose a specific value of given by = 4. For this case, 
the probability density (12. 9p may be explicitly written as 

|Pi,2, -,4( a;i,X2, ■ ■ ■ ,X4 )|^ = exp {2giyi + 2(gi + g2)y2 + 2(gi + g2 + g3)y3} • (2.17) 

Suppose, the conditions (I2.10p for a bound state formation are slightly modified for this 
case as 

gi < 0, gi + g2 = 0, gi + g2 + gs < . 

Taking into account this new condition, it is easy to see that when yi = X2 — xi or 
Vs = X4 — X3 tends towards infinity, the probability density in Eq. (I2.17P still decays 
like a bound state. On the other hand, when y2 = x^ — X2 tends towards infinity, the 
probability density in Eq. (I2.17P remains finite. Hence, the clusters of particles given by 
{1,2} and {3,4} satisfy all the criteria of a clustered state. Generalizing this specific 
example, we replace some of the inequalities in Eq. (I2.10p by equalities. More precisely, 
we modify the conditions in Eq. (I2.10p as 

I 

^^gj = 0, for / G f^AT , 

i=l 

I 

Qi < 0, for / ^{Qn — ^n), 

i=l 

where is any non-empty proper subset of and {fl^ — Q^) is the complementary 
set of Qn- Let us now assume that the quasi-momenta associated with the plane wave 
function (12.70 satisfy the relations (12.80 and (2.18a,b), where the set Qn contains p 
number of elements. Due to Eq. (12.90 it is evident that this plane wave function would 
represent a clustered state containing (p + 1) number of clusters. 

Finally, we try to find the simplest possible condition for which the Bethe state (12.30 
would represent clusters of bound particles. Let us assume that the quasi-momenta 
corresponding to this Bethe state satisfy the relations (12.120 . As a result, the coefficients 
of all plane waves except (12.70 take the zero value within the Bethe state (12.30 . Thus, due 
to the relations (12.120 . the Bethe state (12.30 would reduce to the plane wave function 
(12.70 . Next, we assume that the quasi-momenta corresponding to this Bethe state 



(2.18a) 
(2.186) 
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also satisfy the relations fl2.8p and (2.18a,b). Hence, the plane wave function (12.71) 
represents a clustered state. Consequently, Eqs. (12.81) . fl2.12p and (2.18a,b) together 
yield a sufficient condition for which the Bethe state fl2.3p would represent clusters of 
bound particles. It is obvious that such a Bethe state is expressed through a single 
clustered state, instead of a linear superposition of several clustered and bound states. 

We would like to make a remark at this point. It may be noted that, for the 
case of the (5-function Bose gas, it is not possible to find any set of quasi-momenta 
which simultaneously satisfy the equations fl2.8p . fl2.12p and (2.18a,b). Indeed we have 
already mentioned that, for this case, the conditions fl2.8p and (I2.12p completely fix 
the corresponding quasi-momenta as equidistant and refiection symmetric points on a 
straight line parallel to the imaginary axis in the complex momentum plane. It is easy 
to check that such quasi-momenta do not satisfy Eq. (2.18a). To bypass this problem, 
it should be noted that Eqs. (12. 8p . (I2.12p and (2.18a,b) together yield only a sufficient 
condition, but not the necessary condition for which a Bethe state leads to clusters of 
bound particles. In fact, it is possible to take a partially relaxed version of Eq. (12.120 
given by 

A{kr,kr+i) = 0, for r e n'j^ , (2.19) 

where Q'j^ is any non-empty proper subset of Q^. However, the corresponding Bethe 
state (12. 3 p not only contains the plane wave (12. 7p . but also plane waves like (12. 4p with 
uj taking values within a subset of all possible nontrivial permutations. The Bethe state 
(12. 3p would represent clusters of bound particles if each of these plane waves associ- 
ated with nontrivial permutations behaves like either a bound state or a clustered state. 
Hence, as a consequence of taking Eq. (12.190 instead of Eq. (I2.12p . Eq. (12. lip or ana- 
logues of Eqs. (2.18a,b) for the required values of nontrivial permutation u must also 
be satisfied. For the case of the 5-function Bose gas, a solution of Eqs. (12. 8 p and (I2.19P 
yields N number of quasi-momenta which can be represented through refiection sym- 
metric points on several (more than one) 'strings' parallel to the imaginary axis in the 
complex momentum plane. It can be shown that, apart from satisfying Eqs. (2.18a,b), 
those quasi-momenta also satisfy Eq. (12. lip or analogues of Eqs. (2.18a,b) for the re- 
quired values of nontrivial permutation u. As a result, the Bethe state corresponding to 
those quasi-momenta represents clusters of bound particles for the case of the (5-function 
Bose gas. Evidently, such a Bethe state is expressed through a superposition of several 
clustered and bound states. 

In the next section our aim will be to show that, unlike the case of the (5-function 
Bose gas, it is possible to simultaneously solve the equations (12. 8p . (I2.12p and (2.18a,b) 
for the case of the derivative 5-function Bose gas. To this end, let us discuss how to 
simplify these three equations for the latter case. Since the exact form of the set Qn 
appearing in (2.18a,b) may depend on the value of (p for the case of the derivative 5- 
function Bose gas, in the following we shall replace the notation fl]y hj ^N,(f>- Using 
Eqs. (12. 8p and (I2.12p along with the form of A{ki, km) given in (12. 5p it is easy to see 
that, similar to the case of a localized bound state, the quasi-momenta associated with 
clusters of bound particles can be written in the form (12.130 and the imaginary parts of 
these quasi-momenta satisfy the relation (I2.15p . By using Eq. (I2.15p . we can recast the 
remaining conditions (2.18a,b) for the formation of clusters of bound particles (for any 
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given values of (j) and A^) as 



X ^^^^ sin[(Ar -/)</)] = 0, for I e ^ , (2.20a) 
sm 

^sm{l^ > 0, for / e{QN-^N^), (2.206) 

sm (p 

where x is any non-zero real number, QN,(f> is any non-empty proper subset of and 
{Qn — ^Af,</)) is the complementary set of flN,(i,- Let us assume that Eqs. (2.20a,b) are 
satisfied for some values of and A^, where flN,<f> is given by 

^N,(f> = {h, h, - ■ ■ , Ip} , (2.21) 

with 1 < p < — 1. Then from Eq. f l2.9p it follows that the sets of particles given by 

{1, • • • , h}, {/i + 1, ■ ■ ■ , k}, , {Ip-i + 1, ■ ■ • , Ip}, {Ip + !,■■■ ,N} represent {p + 1) 

number of clusters of bound particles. Moreover, the numbers of particles present within 
each of these clusters, i.e., the size of the clusters, may be written in the form 

{{h,l2-h,--- Jp- Ip^u N-lp}}. (2.22) 

Since the quasi-momenta associated with localized both bound states and clusters 
of bound particles are given by Eq. f l2.13p . the momentum and energy eigenvalues for 
clusters of bound particles can be derived in exactly the same way as has been done 
earlier [A2\ for the case of a bound state. Inserting the quasi-momenta given in Eq. fl2.13p 
to Eqs. (2.6a,b), we obtain the momentum eigenvalue as 

P ^ fix . (2.23) 

sm (p 



and the energy eigenvalue as 



E = — — A— — ^ . 2.24 

sin 20 ^ ' 



3 Farey sequences and clusters of bound particles 

In this section, we shall try to find all possible solutions of the sufficient conditions 
(2.20a,b) for constructing clusters of bound particles in the case of the derivative 6- 
function Bose gas. Some properties of the Farey sequences jl5] in number theory will 
play a crucial role in our analysis. Due to the existence of the parity transformation, as 
mentioned in the earlier section, it is sufficient to concentrate on values of lying in the 
range < < |. Within this range of 0, sin0 > and hence the conditions (2.20a,b) 
for forming clusters of bound particles reduce to 

X sin(/0) sin[(A^ - 00] = , for I e Qn,,p , (3.1a) 
X sin(/0) sin[(A^-/)0] > 0, for / e (n^ - nN,<p) ■ (3.16) 
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Let us first try to find the values of (j) for wliicli Eq. (3.1a) liolds true. For any given 0, 
tliis equation would be satisfied for some value of I if either sin /0 = or sin[(A^— /)</>] = 0. 
For sin/0 = 0, /0 = fevr and so (^/tt = k/l where k is an integer. Since / G ^N,ij), the 
denominator in the expression of (p/rc (= k/l) is always less than A^. Similarly, for 
sin[(A^ — = 0, (p/n = m/{N — I) where m is an integer. Since / G flN,(p, in this 
case also, the denominator in the expression of (p/ir (= m/{N — I)) is less than N for 
all values of /. Hence it follows that, the condition (3.1a) would be satisfied if and only 
if (p/iT can be expressed in the form 

vr 

where {a, 6} are relatively prime integers (i.e, the greatest common divisor of a and h is 
1), taking values within the ranges 

0<a<^, 2<6<A^-1. (3.3a, 6) 

Due to Eq. (3.3b) it is evident that, clusters of bound particles can exist only for > 4. 
In the following, we shall establish a connection of the fractions 0/7r, given by Eqs. (13. 2p 
and (3.3a,b), with the Farey sequences in number theory. 

For a positive integer iV, the Farey sequence is defined to be the set of all the fractions 
a/b in increasing order such that (i) < a < 6 < iV, and (ii) {a, h] are relatively prime 
integers [45] . The Farey sequences for the first few values of are given by 



'2 • 



'5 • 



1 

T T 
oil 

12 1 

112 1 

1 3 2 3 1 

1112 3 1 

1 4 3 2 3 4 1 

1 1 1 2 1 3 2 3 4 1 
----------- (3.4) 

15435253451 ^ ' 



These sequences enjoy several properties, of which we list the relevant ones below. 

(i) Let a/h,a' /h' are two fractions appearing in the Farey sequence F/y. Then a/h < 
a' /b' ( a' /b' < a/b ) are two successive fractions in Fjy, if and only if the following two 
conditions are satisfied: 

a'b - ab' = 1 (-1) , (3.5a) 
b + b' > N. (3.56) 

It then follows that both a and b' are relatively prime to a' and b. 

(ii) For A^>2, ifn/A^isa fraction appearing somewhere in the sequence F/v (this 
implies that and n are relatively prime according to the definition of -F/v), then the 



11 



fractions ai/hi and 02/62 appearing immediately to the left and to the right respectively 
of n/N satisfy 



Oi , a2 < n , and ai + 02 = n , 

bi , b2 < N , and bi + b2 = N . (3.6) 

To apply the above mentioned Farey sequence in the present context, let us define a 
subset of F/v as 



a ^ a 1 , , 



and a subset of FjJ^r as 



^ e n} . (3.* 



Using these definitions of various subsets of a Farey sequence, we find that 

F'j, = F'^_, U . (3.9) 

Furthermore, it is worth noting that, Eqs. (13.21) and (3.3a,b) can equivalently be ex- 
pressed as 

^ G . (3.10) 

TT 

Consequently, it follows that the condition (3.1a) for cluster formation is obeyed if and 
only if ^/tt G F'j^^i- In this context it may be observed that, due to Eq. (13.91) . all the 
elements of -Fjv-i ^Iso present in F'^. By using such an embedding of F'^_^ into F^, 
we find that any fraction a/b & F'n-i belongs to one of the four distinct classes, which 
are defined in the following: 

I. At least one of the fractions nearest to a/6 (from either the left or the right side) in 
the sequence F'^ lies in the set F'^. Then, from a property of the Farey sequences, it 
follows that {6, A^} are relatively prime integers in this case. 

II. None of the nearest fractions of a/b (from the left or right side) in the sequence F'^ 
lies in the set F^, and {6, iV} are relatively prime integers. 

III. N is divisible by b. Clearly, {6, iV} are not relatively prime integers in this case. 

IV. N is not divisible by 6, and {6, A^} are not relatively prime integers. 

To demonstrate the above mentioned classification through an example, let us choose 
A^ = 6. For this case, the sets F5, Fg and Fg' are given by 



1112 
5 4 3 5 



F'. Ill 12. 1 
^ ■ 6 5 4 3 5' ^ ■ 6 



Using the embedding of Fg into Fg, it is easy to verify that each fraction in F5 falls 
under one of the four classes discussed above. More precisely, the fractions 1/5, 1/4, 
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1/3 and 2/5 belong to type I, type IV, type III and type II respectively. Returning 
back to the general case we note that, for any fraction a/h G -Fjv-i ' either 
relatively prime integers or not relatively prime integers. If {6, A^} are relatively prime 
integers, then it is obvious that a/h must be an element of either type I or type II. On 
the other hand, if {6, N} are not relatively prime integers, then a/h must be an element 
of either type III or type IV. In this way, one can show that any fraction a/6 G -F^v-i 
belongs to one of these four distinct classes. This type of classification for the elements 
of 

F'n-i "^ill shortly play an important role in our analysis on the formation of clusters 
of bound particles. 

Next, we try to find the elements of -F^r_i which would satisfy the remaining condition 
(3.1b) for cluster formation. For any a/h G Fjy_i, one can express the rational number 

Na/h as 



Na 



Na 



(3.11) 



where [x] denotes the integer part of x and p G {0, 1, ■ ■ • 6 — 1}. The above equation can 
also be written in the form 

aN-ht = 5, (3.12) 
where the integers t and 5 are defined as 

^Na 



h 
Na 



+ 1, S = p-h, 



■f ^1. 

6 - 2 ' 
if ^>i. 



(3.13) 



One can derive several bounds on the values of t and 6. By using Eq. fl3.13p . it is easy 
to show that 

1*1 - ' (3.14) 



< 

h - 2 



Using Eqs. (3.3a,b) we obtain Na/h > N/{N — 1) > 1, which in turn yields t > 1. 
Moreover, with the help of Eqs. fl3.12p and fl3.13p . we also find that 



t 1 
iV ^ 2 ■ 



(3.15) 



The derivation of Eq. (13.1 5p is given in Appendix A. 

It may be observed that, for any given a/h G -Fjv-i' ^^e values of t and 5 are uniquely 
determined through Eqs. (13. lip and (I3.13p . Thus, we get a mapping of the form 



(a,h,N)^it,5). 



(3.16) 



In the following, we shall show that this mapping considerably simplifies the analysis 
of Eq. (3.1b) by casting it in an alternative form. To this end, let us define a function 
/(/,Ar,0)as 



/(/, AT, 0) = X sin /0 sin(A^ - /)0 . 



(3.17) 
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Substituting the expression of (p given in Eq. (13. 2p to the above equation, we get 

/(/,iV,0)=xsin(^) sin{^(iV-/)}. 
Next, by using the relation fl3.12p . we express /(Z, A^, cj)) in the form 

/(/, N. ^) = (3.18) 
Let us now define a function g{x, 0) as 



b J \ b \ a 



gix,<P)=sm(^-^) , (3.19) 
where x denotes a continuous real variable. Then Eq. (I3.18P can be written as 

/(/, N, 0) = xi-iy^'gil, <P)9il - A, 0) , (3.20) 
where the parameter A is given by 

A = -. (3.21) 
a 

Note that the zero values of the function g{xj(f)) in the variable x are given by 

TTlh 

x = — , m = 0,±l,±2,--- . (3.22) 
a 



Let us consider a value of / such that /(/, A^, 0) 7^ 0, i.e, / G Qj\f — Q]\f^^. Due to Eq. (I3.20p 
it is evident that both g{l, 0) and g{l — A, 0) take nonzero values for such /. Since g{x 
is a continuous function of the variable x, we can write 



,0) 



sgn [g{l - A, 0)] = {-lY^'hgn [g{l, 0)] , (3.23) 

where sgn denotes the sign function and p{l) represents the number of zero points of the 
function g{x, 0) within the interval I — A<x<l{l<x<l — A) when A > (A < 0). 
It is evident that p(/) = for all values of I G Qn — ^N,(t>, when A = 0. Now, by using 
Eqs. fl3:20|l and ([323]), we find that 

sgn [/(/, N, 0)] = sgn [x(-l)*+^+''«] . (3.24) 

With the help of the above equation, we can express Eq. (3.1b) in an alternative form 
given by 



sgn 



for all I enN- n^,^ . (3.25) 



Due to Eq. fl3.22p . it follows that the distance between two consecutive zero points of 
the function g{x, 0) is given by b/a. Moreover, by using Eqs. (I3.14p and (I3.2ip . we find 
that 

- > 2| A I . (3.26) 
a 
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Hence, the value of p{l) can be either or 1 for each /. Since the parameter x does not 
depend on the value of /, it is evident that Eq. fl3.25p would be satisfied if either 



p(/) = , for all IeQn- ^n,^ , (3.27) 

and X is chosen such that sgn{x) = 

p{l) = 1 , for all leQN- ^N,<t> , (3.28) 

and X is chosen such that sgn{x) = (— !)*• However, for any a/b G -Fjv-d i^ '^^^ ^e 
shown that there exists at least one / G — f^^v,,/, for which p{l) = 0. The proof of this 
statement is given in Appendix B. Thus, it is never possible to satisfy the condition given 
in Eq. (13.281) . Consequently, for any given a/b G F^__-^, clusters of bound particles can be 
obtained if and only if Eq. (13.271) is satisfied and x is chosen such that sgn{x) = (~1)*^^- 
On the other hand, clusters of bound particles cannot be formed if Eq. fl3.27p is violated, 
i.e., if the following condition is satisfied: 



p(/) = 1 , for at least one / G VL^ — ^n,iI> ■ (3.29) 

Previously, we have shown that all elements of F'j^_^ can be divided into four distinct 
classes. In the following, we shall analyze each class separately and examine whether 
the elements belonging to each class satisfy Eq. fl3.27p or Eq. fl3.29p . 



Analysis of Case I: Let (p/ir = a/b be a fraction of type I within the set -Fjv-i- 
this case, there exists a fraction n/N & Fj^, such that 



n a 
N b 



or 



a n 
b N 



Hence, using the property of Farey sequences given in Eq. (3.5a), we get 

aN -bn = ±l. 
Comparing the above equation with f l3.12p . we find that 

1 



6 
b 



n-t± 



Combining Eqs. (I3.3ip and (I3.14p . one obtains an inequality of the form 



n-t± 



< 



(3.30) 



(3.31) 



(3.32) 



Since, due to Eq. (3.3b) it follows that 1/6 < 1/3, one can easily show that Eq. ( 13.32P 
would be satisfied if and only if t = n. Since {n, N} are relatively prime integers, it 
follows that {t, N} are relatively prime integers for all fractions of type I. Substituting 
n in the place of t in Eq. fl3.12p . we get 



aN — bn = 6 . 



(3.33) 
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Comparing Eqs. ( I3.30p and fl3.33p . we find that for any fraction of type I, the value of 
S is given by 



±1 



(3.34) 



Next, for the sake of convenience, we consider only those fractions of type I, which 
yield 5 = +1. By using Eq. f l3.2ip . we obtain the corresponding value of A as 



A = -. 

a 



(3.35) 



Let us take any value of I such that /(/, A^, 0) ^ 0, i.e., / G fijv — fi^r,,/,. Using Eqs. fl3.20p . 
(I3.22P and (13.351) . we find that this / satisfies the relations 



mb 
a 



1 mb 
a a 



(3.36) 



for any integer value of m. In the following, we shall try to find if any zero points of the 
function g{x, (p) exists within the range I — 1/a < x < I, i.e., whether mb/a G (Z — 1/a, Z) 
for any integer value of m. To this end, we express mb/a in the form 



mb 
a 



mb 



+ 



(3.37) 



where c G {0, 1, 2, ■ ■ ■ , a — 1}. Let us now assume that mb/a G {I — 1/a, I). Then, by 
using Eq. fl3.37p . we get the relation 



/ < 



mb 



< I 



c 
a 



(3.38) 



However, for c G {0, 1, 2, ■ ■ ■ , a — 1}, one also finds that 

/_1</_1±^, /--</. 
a a 



Combining the above inequalities with those given in Eq. f l3.38p . we obtain the relation 

mb 



l-l< 



<l 



which evidently leads to a contradiction. Hence it is established that, for any fraction of 
type I which gives 5 = 1 and for any / G i^Tv^^Af,^, the condition mb/a G (/ — 1/a, /) can 
never be satisfied. As a result, Eq. f l3.27p is obeyed for all fractions of type I which yield 
5 = 1. Similarly, it can be shown that Eq. fl3.27p is also obeyed for all fractions of type 
I which yield 5 = —1. Consequently, we find that clusters of bound particles are formed 
for all fractions of type I. We have already seen that the relation t = n holds for all of 
these fractions. Therefore, according to the discussion just below Eq. fl3.27p . the sign of 
X should be chosen for this type of clusters of bound particles as sgn{x) = (—1)"^^. 
Analysis of Case II: Let us now consider the fractions of type II within the set -Fjv-i- 
At first our aim is to show that 



|5|>1 



(3.39) 
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for all fractions of this type. By using Eqs. (13. lip and fl3.13p . one can calculate the 
values of t and 6 for these fractions. In contrast to the case of fractions of type I for 
which {t, N} are always relatively prime integers, for the present case {t, N} may or 
may not be relatively prime integers. As an example, let us take N = 7, for which 
a/b = 1/5 is a type II fraction. Since 7/5 = 1 + 2/5, using fl3.13p we find that t = 1 and 
6 = 2. Hence, {t, N} are relatively prime integers in this case. On the other hand, we 
may also take = 6, for which a/b = 2/5 is a type II fraction. Since 6x2/5 = 2 + 2/5, 
it follows that t = 2 and 6 = 2. Hence, {t, N} are not relatively prime integers in this 
case. 

Let us first consider all fractions of type II within the set -F/v-i? which {t,N} 
are relatively prime integers. As we have seen that t > 1 and Eq. f l3.15p also gives 
t/N < 1/2, it is evident that 

^ e f;^ c n . (3.40) 

However, since a/b is a fraction of type II, a/b and t/N cannot be two consecutive 
fractions in the sequence F'j^. Therefore, at least one of the two relations given in 
Eqs. (3.5a) and (3.5b) must be violated if we choose a' = t, b' = N. Since Eq. (3.5b) is 
obviously satisfied for this case, Eq. (3.5a) must be violated. As a result, we obtain 

\aN-bt\j^l. (3.41) 

Next, let us assume that aN — bt = 0. Since {a, b} and {t, N} are both relatively prime 
pairs, the relation aN — bt = implies that a = t and b = N. However, since a/b G -Pjv-i; 
the relation b = N evidently leads to a contradiction. Thus, it follows that 

aN -btj^O. (3.42) 

Since a, b, t, N are all integer numbers, combining Eqs. fl3.4ip and fl3.42p we find that 

\aN -bt\>l. (3.43) 

Using Eqs. fl3.43p and fl3.12p . one can easily establish the validity of Eq. (13.390 . 

Next, we consider all fractions of type II within the set -F/v-i? which {t,N} are 
not relatively prime integers. In this case, one can express t and A^ as 

t = at' , N = aN' , (3.44a, b) 

where {t', N'} are relatively prime integers and a (> 1) is an integer. Using Eqs. (13.120 
and (3.44a,b), we obtain 

\5\ = a\aN' -bt'\. (3.45) 

Let us now assume that aN' — bt' = 0, i.e, a/b = t'/N'. Since {a,b} and {t',N'} are 
both relatively prime pairs, the relation a/b = t'/N' implies that a = t' and b = N' . 
Thus, by using Eq. (3.44b) we obtain A^ = ab. However, the relation N = ab clearly 
contradicts our initial assumption of a/b being a fraction of type II, for which A^} 
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must be relatively prime integers. Consequently, we find that aN' — ht' ^ 0. Since 
a, 6, t', A^' are all integer numbers, it follows that 



\aN' -ht'\ > 1. (3.46) 

Using Eqs. f l3.45p . fl3.46p and the relation a > 1, one can easily establish the inequality 
given in f l3.39p . 

Next, our aim is to show how Eq. (13.291) follows from the inequality (13.391) . which 
is obeyed by all fractions of type II. To this end, let us assume that c/d is the nearest 
fraction of a/b in the sequence F/v, either from the left side, i.e., 

Fn: , ^, ^, (3.47a) 



or, from the right side, i.e., 

Fn: ' ^ ^ ' ^3.476) 

Since a/b is a fraction of type II, it is evident that d < N —1 (note that this statement 
remains valid for any a/b which is not a fraction of type I). By using the method of 
contradiction, in the following we shall show that d G Qn — ^n,(I> for both of the cases 
(3.47a) and (3.47b). Let us first assume that a; = ci is a zero point of the function 
g{x, (p). Thus, by using (I3.22p we obtain d = mb/a, i.e, 

m = —d . 
b 

Since {a, b} are relatively prime integers, the above equation can only be satisfied if 
d = ab, where a is a positive integer. However the relation d = ab leads to a contradic- 
tion, since the denominators of two consecutive fractions in a Farey sequence must be 
relatively prime integers. Hence, we find that x = d is not a zero point of the function 
g{x, 0). Next we assume that x = d is a zero point of the function g{x — A, 0). Thus, 
by using (I3.22p we obtain d — 6/a = m'b/a, i.e., 

ad — m'b = 6 . 

Equating the (left hand side) l.h.s. of the above equation with that of Eq. (I3.12p . we 
get 

a t — m' 



b N-d 

Since {a, b} are relatively prime numbers, the above relation implies that t — m' = Pa 
and N — d = (3b, where /3 is a positive integer. However, it is easy to see that the 
relation N — d = Pb contradicts Eq. (3.5b), which is satisfied by the denominators of 
two consecutive fractions in a Farey sequence. Hence, we find that x = d is not a zero 
point of the function g{x — A, 0). Since d < N — 1 and x = d is not a zero point of 
either g{x, 0) or g{x — A, 0), it is established that d ^ Qn — ^N,(f> for both of the cases 
(3.47a) and (3.47b). 
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Due to Eq. fl3.39p . all fractions of type II can be subdivided into two classes charac- 
terized by 5 > 1 and 6 < —1 respectively. Let us first consider all fractions of type II 
for which 6 > 1. For this case, we assume that c/d is the nearest fraction of a/b (in the 
sequence F^) from its left side, as shown in Eq. (3.47a). Consequently, by using (3.5a), 
we obtain 

ad-bc = l. (3.48) 

By choosing m = c in Eq. 03.221) . one finds that x = c6/a is a zero point of the function 
g{x,(l)). By using Eq. 03.481) . we calculate the difference between d and cb/a as 

/ , cb\ , c6 1 , , 

Did, — =d = -. (3.49) 

\ a J a a 

Since we are considering the case 6 > 1, Eq. f l3.2ip yields A > 1/a. As a result, 
Eq. 03.491) leads to an inequality given by 

D(^d, <A, (3.50) 

which shows that the zero point x = cb/a of the function g{x, 0) lies within the interval 
{d — A,d). Consequently, it is established that p{d) = 1 for all fractions of type II with 
^ > 1. 

Next, let us consider all fractions of type II for which 6 < —1. For this case, we 
assume that c/d is the nearest fraction of a/b (in the sequence F^) from its right side, 
as shown in Eq. (3.47b). Consequently, by using (3.5a), we obtain 

ad — be = —1 . (3.51) 

By using the above relation, we obtain the distance between cb/a and d as 

d('1, d]='^-d='i. (3.52) 
\ a J a a 

Since we are considering the case 6 < —1, Eq. 03.2ip yields l/a<— A. As a result, 
Eq. 03.521) leads to an inequality given by 

D (^^, d^ <-A, (3.53) 

showing that the zero point x = cb/a of the function g{x,(f)) lies within the interval 
{d, d — A). Consequently, it is established that p{d) = 1 for all fractions of type II with 
S < -1. 

From the above analysis, it is clear that the condition 03.29P is obeyed for all fractions 
of type II. Consequently, clusters of bound particles cannot be formed for any fraction 
of this type. 

Analysis of Case III: Let (p/n = a/b be a fraction of type III within the set -Fjv-i- 
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Since N is divisible by b in this case, we may write = pb, where p (> 1) is an integer. 
Using Eqs. f l3.1ip and f l3.13p . we find that t = ap and 5 = for this case. Substituting 
5 = in Eq. (13.211) . we obtain A = 0. Hence, due to Eq. (13.23 p . it trivially follows 
that p{l) = for all values of I G — ^N,<f>, i-e., Eq. (I3.27P is satisfied for this case. 
Therefore, we find that clusters of bound particles are formed for all fractions of type 
III. According to the discussion just below Eq. (I3.27p . the sign of x should be chosen 
for this type of clusters of bound particles as sgn{x) = (—1)"^"*'^. 

Analysis of Case IV: Let (p/n = a/b be a fraction of type IV within the set i^jy-i- 
For this case, we can write 

b = ab', N = aN', (3.54a, 6) 

where a, b', N' are some integers such that a > 1, b' > 1 and {b', N'} are relatively 
prime integers. Using Eqs. (13.120 and (3.54a,b), we obtain 

S = a{aN' - tb') . (3.55) 

Let us now assume that aN' — tb' = 0, which gives 

a = §j. (3.56) 

Since {b', N'} are relative prime integers, the above equation would be satisfied if t = 
/3N', where /3 is a positive integer. Substituting this value of t in Eq. (I3.56p . we obtain 

a = (3b'. (3.57) 

Eqs. (3.54a) and (I3.57P imply that {a, b} are not relatively prime integers, which contra- 
dicts our basic assumption that a/b G -Fiv-i- Thus it is established that aN' — tb' ^ 0. 
Since a, V, t, 6' are all integers, it follows that 

\aN'-ty\>\. (3.58) 

Combining Eqs. (I3.55P and (I3.58P we find that, similar to the case of fractions of type 
II, the inequality given by 

|5|>1, 

is satisfied for all fractions of type IV. Consequently, by carrying out the rest of the 
analysis in exactly same way as we have done in this section for the fractions of type II, 
it can be shown that all fractions of type IV obey Eq. (13.290 . Hence, clusters of bound 
particles cannot be formed for any fraction of this type. 



4 Some properties of clusters of bound particles 

In the previous section we have shown that, for the case of derivative the 5-function 
Bose gas with given values of and N , the Bethe state (12. 3p represents clusters of 
bound particles if ^/tt is either a fraction of type I or type HI within the set -Fjv-i- 
In this section, our aim is to find the number of clusters present within such a Bethe 
state and the sizes of these clusters (i.e., number of bound particles present in each of 
these clusters). We would also hke to investigate the behavior of these clusters of bound 
particles under small variations of the coupling constant. 
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4.1 Sizes of the clusters of bound particles 

Let us first consider clusters of bound particles when ^/vr = a/h is taken as any fraction 
of type I within the set F'j^_^. In section 2 we have seen that, to find the number and 
sizes of the clusters within a Bethe state, we have to determine the set Qn,<I) for which 
Eq. (3.1a) is satisfied. Since {A^, b} are relatively prime integers for any fraction of type 
I, we can express as 

N = pb + r, (4.1) 

where 1 < r < 6 — 1. Hence, for the discrete variable I taking values within the set Qn, 
the zero points of the functions sin/</) and sin(A^ — /)0 are respectively given by the sets 

Si = {b,2b, ,pb}, S2 = {N -b,N -2b, ,N-pb}. (4.2a, 6) 

Combining the sets Si and 5*2 by using Eqs. (14. ip and (4.2a,b), we obtain flN,<f> as 

Qn,^ = SiUS2 = {r, b, r + b, 2b, , r + {p - 1)6, pb} . (4.3) 

Comparing (14. 3 p with (12. 21 p . and also using (I2.22p . it is easy to see that the sizes of the 
clusters are given by 

{{r, b — r, r, b — r, , r, b — r, r}} . (4.4) 

Hence, for any fraction of type I, the corresponding Bethe state contains {p+ 1) number 
of clusters of size r and p number of clusters of size [b — r). Next, by using the method 
of contradiction, we would like to show that these two possible sizes of the clusters, i.e., 
r and (6 — r), must be relatively prime integers. To this end, let us first assume that b 
and r are not relatively prime integers. Therefore, we can write b and r as 6 = ab' and 
r = ar', where a > 1. Substituting these values of b and r in Eq. (14. ip . we find that 

N = a{pb' + r') . 

Thus a is a common factor of and b. However, this result contradicts the fact that 
{A^, 6} must be relatively prime integers for any fraction of type I. Hence it is established 
that {b, r} are relatively prime integers. From this relation, it trivially follows that 
{r, 6 — r} are relatively prime integers and, in particular, r ^b — r. Consequently from 
Eq. (14. 4p we find that, for any fraction of type I, the corresponding Bethe state contains 
heterogeneous clusters of two different sizes. As a special case, let us consider any 
fraction of type I with denominator satisfying the relation b > N/2. Due to Eq. ( 14. ip it 
follows that, p = 1, r = N — b and b — r = 2b — N for this case. Hence, the corresponding 
Bethe state contains two clusters of the size {N — b) and one cluster of the size {2b — N). 

Next, we consider the clusters of bound particles corresponding to any fraction of 
type HI. In this case A^ can be written as A^ = pb, where p is an integer greater than 
one. Consequently, for the variable / taking value within the set Qm, the zero points of 
the functions sin/0 and sin(A^ — /)0 coincide with each other and yield flN,(t> as 

n^,^ = {b, 2b, 3b, ,{p- 1)6} . (4.5) 
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Hence p number of clusters are formed in this case. Comparing (14 .Sp with f l2.2ip . and 
also using fl2.22p . it is easy to see that each cluster of this type has the size b. In 
other words, the corresponding Bethe state fl2.3p contains N/b number of homogeneous 
clusters, each of which is made of b number of bound particles. 

In Table 1 we show all the fractional values of ^/vr for which clusters of bound 
particles exist within the range of N given by 4 < < 10, the types of these fractions, 
the values of S for fractions of type I, and the sizes of the corresponding clusters using 
the notation of Eq. (1^:^ . 





Value of 0/7r 


Type 


Value of 6 


Size of the clusters 


4 


1/3 


I 


1 


{{1,2,1}} 


5 


1/4 


I 


1 


{{1,3,1}} 


5 


1/3 


I 


-1 


{{2,1,2}} 


6 


1/5 


I 


1 


{{1,4,1}} 


6 


1/3 


III 


— 


{{3,3}} 


7 


1/6 


I 


1 


{{1,5,1}} 


7 


1/4 


I 


-1 


{{3,1,3}} 


7 


1/3 


I 


1 


{{1,2,1,2,1}} 


7 


2/5 


I 


-1 


{{2,3,2}} 


8 


1/7 


I 


1 


{{1,6,1}} 


8 


1/4 


III 




{{4,4}} 


8 


1/3 


I 


-1 


{{2,1,2,1,2}} 


8 


2/5 


I 


1 


{{3,2,3}} 


9 


1/8 


I 


1 


{{1,7,1}} 


9 


1/5 


I 


-1 


{{4,1,4}} 


9 


1/4 


I 


1 


{{1,3,1,3,1}} 


9 


1/3 


III 




{{3,3,3}} 


9 


3/7 


I 


-1 


{{2,5,2}} 


10 


1/9 


I 


1 


{{1,8,1}} 


10 


1/5 


III 




{{5,5}} 


10 


2/7 


I 


-1 


{{3,4,3}} 


10 


1/3 


I 


1 


{{1,2,1,2,1,2,1}} 


10 


2/5 


III 




{{5,5}} 



Table 1: The fractional values of ^/vr for which clusters of bound particles exist for 
4 < < 10, the types of these fractions, the values of 6 for fractions of type I, and the 
sizes of the corresponding clusters are shown. 



4.2 Stability of clusters of bound particles and binding energies 

We want to explore here how the nature of clusters of bound particles change under slight 
variations of the coupling constant. For any given value of A^ (> 4), we first choose a 
value of the coupling constant within the range < < | such that clusters of bound 
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particles are formed. If one increases or decreases this value of (f) by an infinitesimal 
amount, it is obvious that all inequalities in Eq. (3.1b) continue to be satisfied, but all 
equalities in Eq. (3.1a) are transformed into some inequalities. Therefore, clusters of 
bound particles cease to exist even for a very small change of the coupling constant. 
One of the following two different cases can occur in such a situation. In the first case, 
at least one of the equalities in Eq. (3.1a) is transformed into an inequality of the form 

X sin(/0) sin[(A^ - /)0] < . (4.6) 

It is evident that, the probability density of the corresponding Bethe state f l2.3p would 
diverge if at least one of particle coordinates tends towards infinity (keeping the centre of 
mass coordinate fixed). As a result, this Bethe state becomes ill-defined and disappears 
from the Hilbert space of the Hamiltonian (11. 2p of derivative the 5-function Bose gas. 
So we may say that clusters of bound particles become unstable in this case. Let us 
now consider the second case, for which all of equalities in Eq. (3.1a) are transformed 
into inequalities of the form 

X sin(/0) sm[{N-l)(f)] > 0, (4.7) 

due to an infinitesimal change of the coupling constant 0. It is evident that, for this case, 
Eqs. (3.1a,b) are transformed to Eq. (I2.16P within the range of given by < < ^. 
As a result, clusters of bound particles merge with each other and produce a localized 
bound state containing only one cluster of particles. Therefore, in this case, we may say 
that clusters of bound particles turn into a localized bound state with only one cluster. 

In this context, we note that the first and second derivatives of the function f{l,N, 0) 
in Eq. (I3.17P are given by 

— = — sm(iV0) ^ sm(A^ - 2/)0 , (4.8a) 

— = cos{N(j)) cos{N - 2/)0 . (4.86) 

By calculating these derivatives at some value of for which clusters of bound states 
exist and at a value of / within the set ^lN,<j>, one can find whether the corresponding 
equality in Eq. (3.1a) is transformed into an inequality of type (14. 6 p or (14. 7p due to 
an infinitesimal change in the coupling constant. In the following, we shall use this 
procedure to determine the nature of clusters of bound particles for any given values of 
and N. 

Let us first consider clusters of bound particles when (p/n = a/b is taken as any 
fraction of type III. As mentioned earlier, in this case N can be written as N = pb (where 
p is an integer) and from Eq. (I3.13P it follows that t = ap. Moreover, due to Eq. (14. 5p . 
/ should be taken to be of the form I = mb, where m is an integer. For the values of N, 
and I given by iV = pb, = ira/b and I = mb, we get sin(A^0) = sin(7rap) = and 
sin(A^ — 2/)0 = sin(A^0 — 27rma) = 0. Thus, by using Eq. (4.8a), we obtain 
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which shows that (p = ira/b is an extremum point of the function /(/, A^, (p). To determine 
the nature of this extremum point, it is necessary to calculate the second derivative of 
f{l, N, 0). To this end we recall that, according to our discussion just after Eq. fl3.27p . 
the value of x should be chosen such that the inequality (— < is satisfied. Since 
t = ap for all fractions of type III, we can rewrite the above mentioned inequality as 

i-irx < 0. (4.10) 

Moreover, for the case = pb, = 7ca/b and / = mb, we obtain cos(A^0) = cos(A^ — 
2/)0 = (—1)"^. Substituting these values of cos(A^0) and cos(A^ — 2/)0 in Eq. (4.8b), 
and subsequently using the inequality (14.101) . we find the relation 



2liN - I) . xi-ir < , (4.11) 



O ' 



which shows that the function f[l,N,(j)) has a maxima at = na/b. Consequently, if 
one moves slightly away from the point (j) = ira/b {in positive or negative direction), the 
corresponding equality in Eq. (3.1a) is transformed into an inequality of the form fl4.6p . 
Therefore, we find that the clusters of bound particles are unstable around the point 
= na/b, where a/b is a fraction of type III. 

Next, we consider clusters of bound particles corresponding to any fraction of type 
I. Due to Eq. fl3.34p . it is possible to divide all fractions of type I into two subclasses, 
corresponding to the value of S given by +1 and —1 respectively. We shall see below 
that clusters of bound particles related to these two subclasses behave differently under 
the variation of the coupling constant. For any fraction of type I, the union of the 
sets Si and 5*2 given in Eqs. (4.2a,b) yields the set flN,(t>- Therefore, the values of the 
corresponding I should be taken either as I = mb or a.s I = N — mb, where m is an 
integer. Let us first consider the case when / = mb. Since sin(A^ — 2Z)0 = sin(A^0) in 
this case, by using Eq. fl3.12p we obtain 



sin(iV - 2l)(j) = sm{N(j)) = {-ly sin j . (4.12) 

Similarly, for the case I = N — mb, we obtain 

- sin(A^ - 2/)0 = sin(A^0) = (-1)* sin (^^^ . (4.13) 

Substituting the values of sin(A^0) and sin(A^ — 2/)0 given in Eqs. fl4.12p and fl4.13p 
to Eq. (4.8a) for the cases / = mb and I = N — mb respectively, and also using the 
inequality (— l)*x < 0, we get 



df 



™..X(-I)'.si„(^).{;°[:^:i- (4.14) 



Hence we find that, for the case of fractions of type I with 6 = 1, clusters of bound 
particles would become unstable due to a slight increase of the value of (p. On the 
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other hand, such clusters of bound particles would transmute to a localized bound state 
containing only one cluster if the value of (f) is slightly decreased. For the fractions of 
type I with 5 = — 1, clusters of bound particles would behave in exactly reversed order 
if one slightly increases or decreases the value of (f). 

It is interesting to calculate the binding energies of the clustered states that we 
have been considering so far. The binding energy Eb of a clustered state is defined as 
follows. The momentum P and energy E of such a state were presented in Eqs. (12.231) 
and fl2.24p . Let us now consider a state of N free particles which are all moving with 
the same momentum p so that the total momentum Np is equal to the expression in 
Eq. fl2:23|) . Clearly, 

fc sin(Ar0) 
JS sm0 

The energy of this free particle state is given by 

^0 = Np' = \ . ^ . 4.16 

M sm (0) 

The binding energy Eb corresponding to clusters of bound particles is defined as the 
difference of the energies of the free particle state and the clustered state, i.e.. 



h'x sai{N(j)) 
A^sin^(0) cos( 

h{N,(j)) = X [sin(A^0)cos0 - A^cos(A^0)sin0]. (4.176) 



Eb ^ Eo- E = KNA) , (4.17a) 



It should be noted that Eqs. (4.17a,b) also give the binding energies of localized bound 
states for the case of the derivative 5-function Bose gas [l2llll]. We now consider the 
two kinds of clusters of bound particles in which (p/ir = a/b corresponds to fractions of 
types I and III within the sequence i^iv-i- -^^^ ^^^^ fractions of type III, iV is a 
multiple of b; hence sin(A^0) = 0. Thus the momentum p and the energies E and Eq 
are all equal to zero in this case, and Eb = also. We now turn to fractions of type 
I. In Refs. [121111], it was shown that h{N,(j)) > for all these values of (j). Eq. (4.17a) 
then implies that Eb has the same sign as xsin(A^0). Combining Eqs. (14.12114.1^ . we 
see that x sm{N(f)) and hence Eb are negative for 6=1 (i.e., if {aN — 1) /b is an integer) 
and are positive for 6 = —1 (i.e., if {aN + l)/b is an integer). 

In our earlier works it was shown that, for any given value of > 4, the derivative 
^-function Bose gas allows bound states in only certain non-overlapping ranges of the 
coupling constant called 'bands' (the union of which yields a proper subset of the 
full range of 0), and the location of these bands can be determined exactly |12|I1^. 
Applying the results obtained in this section, we shall now present an alternative method 
of computing the location of these bands within the full interval (0,7r/2) for positive 
values of 0. Since fractions of type I only lead to clusters of bound particles which 
transmute to bound states with a single cluster through infinitesimal changes of in 
appropriate directions, it is evident that any fraction of type I or or 1/2 must appear 
at the two edges of any continuous range of (p/n within which localized bound states are 
formed. Therefore, we include the numbers and 1/2 along with the fractions of type 
I (corresponding to the set -Fjv-i) '^^^^ ^^i^ new set as type I'. Next, we denote two 
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fractions of type I' as ai/hi and a^/h[, and define them as 'conjugate pair' if they appear 
within the Farey sequence F/v as 



Fn- 




(4.18) 



where rii/N & Fj^. Let us also assume that all possible conjugate pairs of this type can 
be obtained by varying the index i within the range i G {1, 2, ■ ■ ■ , a{N)}. Clearly, for 
any > 4, conjugate pairs defined in Eq. (14.181) can be related to three different cases: 
(i) Qi/bi and a[/b[ are both fractions of type I, (ii) ai/bi = and ayb[ is a fraction of type 
I, and (iii) a'Jb'- = 1/2 and Oj/fej is a fraction of type I. Let us first consider any conjugate 
pair corresponding to the case (i). Using Eqs. (3.5a) and (I3.12p along with the relation 
t = n, it is easy to see that S = —1 for ai/bi and 6 = 1 for a'Jb[. Therefore, from the 
discussion of the previous paragraph it follows that clusters of bound particles would be 
transformed to a localized bound state only if one increases the value of (p/ it from a^/bi 
or decreases the value oi (p/ir from a'Jb[ by an infinitesimal amount. Moreover it should 
be noted that, since ai/bi and a'i/b^ are consecutive fractions within the set -F/v-i) ^^e 
sign of the function f{l,N,(j)) cannot flip even for any finite amount of change of (p/n 
within the interval {ai/bi, a'Jb'j). Consequently, for any conjugate pair corresponding to 
case (i), localized bound states continue to exist if ^/vr takes any value within the range 
{ai/bi,a'Jb'^. By repeating similar arguments, one can reach the same conclusion for 
any conjugate pair associated with cases (ii) or (iii). Thus we find that, localized bound 
states can be constructed for any value of within the non-overlapping ranges or bands 
given by 



where i G {1,2, ■■ ■ ,a{N)}. Due to Eq. f l4.18p it is evident that there exists an one- 
to-one correspondence between the conjugate pairs formed by fractions of type I' and 
the elements of the set F^. Therefore, the elements of Fj^ can be used to uniquely 
characterize the bands appearing in Eq. (14.191) . Moreover, it is easy to check that 
Eq. (I4.18P provides the only possible way of relating two fractions like a/b and c/d of 
type I', such that localized bound states are formed for any value of (p/n within the range 
{a/b, c/d). Hence, localized bound states can only be constructed within the ranges of 
(j) given in Eq. (I4.19p . 

For an illustration of the idea of bands, and the stability of clusters of bound par- 
ticles under a variation of and their binding energy discussed above, we calculate 
the binding energies of localized bound states within all bands for = 7, 8 and 9 by 
using Eqs. (4.17a,b), and we show these binding energies as a function of (p/ir in Fig. 1. 
Comparing this figure with the appropriate entries in Table 1, we see that the values of 
(P/tt which are fractions of type III do not lie in any bands, while fractions of type I lie 
at the end points of the bands. Furthermore, fractions of type I with 5 = 1 (—1) lie at 
the right (left) end of the bands and have Eb < (> 0) respectively. 




(4.19) 
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-8 I 

0.1 0.2 0.3 0.4 0.5 

Figure 1: The binding energy Eb of the locahzed bound states as a function of (p/n for 
iV = 7, 8 and 9. 

5 Conclusion 

In this article, we have explored how clusters of bound particles can be constructed in 
the simplest possible way for the case of an exactly solvable derivative 5-function Bose 
gas. To this end, we consider a sufficient condition for which Bethe states of the form 
in fl2.3p would lead to clusters of bound particles. This sufficient condition for obtaining 
clusters of bound particles has not attracted much attention earlier, probably because it 
does not yield any solution at all for the case of the 5-function Bose gas. However we find 
that, for the case of the derivative 5-function Bose gas, this sufficient condition can be 
satisfied by taking the quasi-momenta of the corresponding Bethe state to be equidistant 
points on a single circle having its centre at the origin of the complex momentum plane. 
Furthermore, the coupling constant (0) and the total number of particles (A^) of this 
derivative (5-function Bose gas must satisfy the relations in (3.1a,b). For any given 
> 4, it is found that Eq. (3.1a) is satisfied if (p/ir takes any value within the set -F)v_i, 
which is a subset of the Farey sequence F^-i- Then we classify all fractions belonging 
to the set F'j^_^ into four types. It turns out that fractions of types I and III belonging 
to F'j^_^ satisfy the remaining relation (3.1b), while type II and type IV fractions do 
not satisfy (3.1b). Consequently, clusters of bound particles can be constructed for the 
derivative 5-function Bose gas only for special values of (p/ir given by the fractions of 
types I and III within the set -Fjv-i- 

We also computed the sizes of the above mentioned clusters of bound particles, i.e., 
the number of particles present within each of these clusters. We find that any fraction 
of type I within the set F^_^ leads to heterogeneous clusters of bound particles having 
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two different sizes. On the other hand, any fraction of type III within the set -Fjv-i 
leads to homogeneous clusters of bound particles having only one size. Interestingly, 
clusters of bound particles associated with fractions of type I and type III transform 
in rather different ways under a small variation of the coupling constant. For example, 
it is found that clusters of bound particles associated with fractions of type III cease 
to exist for any small change of the coupling constant. On the other hand, clusters of 
bound particles corresponding to fractions of type I turn into localized bound states 
consisting of a single cluster if the value of the coupling constant is slightly increased or 
decreased. Finally, we find that clusters of bound particles associated with fractions of 
type III have zero binding energy, while clusters associated with fractions of type I have 
either positive or negative binding energy. 

Even though in this paper we have analyzed a particular type of sufficient condition 
for constructing clusters of bound particles in the case of the derivative 5-function Bose 
gas, in future it would be interesting to explore other possible ways of constructing clus- 
ters of bound particles for this exactly solvable system. We have previously mentioned 
that, clusters of bound particles are constructed for the case of the (5-function Bose gas 
by taking the corresponding quasi-momenta to be discrete points lying on several par- 
allel straight lines in the complex momentum plane. In analogy with this case, one may 
try to construct clusters of bound particles for the present case by assuming that the 
corresponding quasi-momenta lie on several concentric circles in the complex momentum 
plane. However it should be noted that, for the case of the derivative 5-function Bose 
gas with any given value of the coupling constant, localized bound states can appear 
only for some particular values of the particle number |32lll3]. Hence, for constructing 
clusters of bound particles in this case, the number of quasi-momenta on each concen- 
tric circle must be taken from the subset of such allowed particle numbers. Due to this 
constraint on the numbers of quasi-momenta, it remains a rather challenging problem to 
construct all possible clusters of bound particles in the case of the derivative (5-function 
Bose gas. 
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Appendix A 

In this Appendix, we shall explicitly derive the relation f l3.15p . At first, we shall show 
that if a/6 G -F^v-i satisfies the condition 

all , , , 

then t/N will satisfy the relation fl3.15p . To this end, we rewrite Eq. f l3.12p as 

N b bN ^ ' 

Eqs. (Al) and (A2) imply that 

t 1 1 5 

N - 2~2N~bN' ^ ^ 

For the case 6 > 0, the relation (13.151) directly follows from Eq. (A3). For the case 6 < 0, 
Eq. (Km yields 

b 2 ^ ' 

Combining (A3) and (A4), it can be easily seen that the relation (I3.15P is satisfied for 
S < 0. Thus we are able to prove that the relation (I3.15P follows from the condition 
(Al) for all possible values of 6. 

Next, for the purpose of proving the validity of condition (Al), let us consider the 
cases of even and odd separately. For the case of even values of A^, we may write 
A^ = 2m. Using the properties (3.5a,b) of Farey sequence, we find that 

m — 1 1 

F2m-1 '■ 



2m - 1 2 

i.e., the fraction nearest to 1/2 from the left side of the sequence F2m-i is given by 
(m — l)/(2m — 1). Consequently, any fraction a/b e -p2m-i satisfies the relation 

a m — 1 , , , 

6 - 2m - 1 ^ ' 

Using (A5), we easily find that the condition (Al) with A^ = 2m is satisfied by any 
a/b e -p2m-i- Next, we consider the case of odd values of A^ and write A^ = 2m + 1. 
Using Eqs.(3.5a,b), it is easy to check that there are two consecutive fractions in the 
Farey sequence given by 

m — 1 1 

F2m '■ 



2m - 1 2 

Consequently, any fraction a/b G also satisfies the relation (A5). Using (A5), we 
find that the condition (Al) with A^ = 2m + 1 is satisfied by any a/b G -^2™- Thus it 
is established that the condition (Al) is satisfied by any a/b G -Fjv-i both even and 
odd values of A^. 
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Appendix B 

Here our aim is to show that for any ajh ^ -^iv-i' there exists at least one I E Qn ~^n,<i> 
such that p(/) = 0. To this end, we shall consider the cases 6 > and 6 < separately. 

Case (i) Let us first consider the case S > 0, for which A > 0. Using Eqs. (3.3a) and 
03 -Mp . we obtain the relation 

0<a + 5<b. (Bl) 

By using the method of contradiction, at first we shall show that (b — 1) E — ^n,4>- 
Let us first assume that a; = 6 — 1 is a zero point of the function g{x,(j)). Then by using 
Eq. fl3.22p we obtain b — 1 = mb/a, which yields 

a 

— = a — m . 
b 

However, the above equation leads to a contradiction because its l.h.s. is a proper 
fraction and its right hand side (r.h.s.) is an integer. Thus x = b — 1 cannot be a zero 
point of the function g{x,(f)). Next, we assume that x = 6 — 1 is a zero point of the 
function g{x — A, (p). Then by using Eq. fl3.22p we obtain b — 1 — 6/a = m'b/a , which 
yields 

a + 6 , 
= a — m . 

b 

Again the above equation leads to a contradiction because its r.h.s. is an integer and 
its l.h.s. is a proper fraction due to Eq. (Bl). Thus x = b — 1 cannot be a zero point of 
the function g{x — A,(f)). Since (6 — 1) < (A^ — 1) and x = 6 — 1 is not a zero point of 
either g{x, 0) or g{x — A, 0), it follows that (6 — 1) G Qn — ^Af,</.- 

Next, our aim is to show that p{l) = for the choice I = b — 1. To this end, let us 
assume that x = m"b/a is a zero point of the function g{x,(f)) and m"b/a < {b — 1). 
Since the difference between {b — 1) and m"b/a may be written as 

ofb-l,^] =-l + -{a-m"), (52) 
\ a J a 

it is clear that m" < a. Using Eqs. (3.3a) and (B2), it is easy to see that D{b — 1, m"b/a) 
yields the minimum positive value for the choice m" = a — 1. Furthermore, by using 
Eqs. (Bl) and fl3.2ip . we find that this minimum positive value of D{b — l,m"b/a) 
satisfies the relation 



, , m"b 
D b-h 



— > A. 



=a-l 



From the above equation it follows that, there exists no zero point of the function g{x, (j)) 
in the variable x within the interval (6 — 1 — A, 6 — 1). Thus, for any a/b G i^^v-i with 
5 > 0, we are able to establish that p{b — 1) = 0. 

Case (ii) Let us consider the case 5 < 0, for which A < 0. Using Eqs. (3.3a) and 
fl3.14p . we obtain the relation 

0<a-5<6. (B3) 
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By using the method of contradiction, in the following we shall show that (6 + 1) G 
f^Tv — ^N,<i>, i-e, (6+1) < (A^ — 1) and a; = (6 + 1) is not a zero point of either g{x, (p) 
or g{x — A,(f)). For this purpose it may be noted that , since a/b e -Fjv-i' "^^^ "^^^ 
principle choose 6 = A^ — 1, which contradicts the condition (6+1) < (A^ — 1). However, 
by applying Eqs. fl3.1ip and fl3.13p for this case, we obtain [A^a/6] = a and 6 = a > 0. 
Thus the choice 6 = — 1 is not compatible with the present case associated with 
negative values of S. Next we assume that x = 6 + 1 is a zero point of the function 
g{x, 0). Then by using Eq. f l3.22p we obtain 6 + 1 = mb/a, which yields 

a 

— = m — a . 
6 

However, the above equation leads to a contradiction because its l.h.s. is a proper 
fraction and its r.h.s. is an integer. Thus x = 6 + 1 cannot be a zero point of the function 
g{x,(f)). Next, we assume that x = 6 + 1 is a zero point of the function g{x — A,0). 
Then by using Eq. f l3.22p we obtain 6 + 1 — 5/a = m'b/a, which yields 

a — 6 , 

— ; — = m — a . 



Again, the above equation leads to a contradiction because its r.h.s. is an integer and 
its l.h.s. is a proper fraction due to Eq. (B3). Thus x = 6 + 1 cannot be a zero point of 
the function g{x — A,(j)). As a result, we find that (6+1) G Qn — ^In,<i> for any fraction 
a/b G -Fjv-i which leads to 5 < 0. 

Next, our aim is to show that p{l) = for the choice / = 6 + 1. To this end, let 
us assume that x = m"b/a is a zero point of the function g{x, (p) and m"b/a > (6+1). 
Since the difference between m"b/a and (6 + 1) may be written as 

D(^,b + l] =-l + ^(m"-a), (BA) 



it is evident that m" > a. With the help of Eqs. (3.3a) and (B4), it is easy to check 
that D{m"b/a,b + 1) yields the minimum positive value for the choice m" = a + 1. 
Furthermore, by using Eqs. (B3) and (13.210 . we find that this minimum positive value 
of D{m"b/a ,6+1) satisfies the relation 



D 



( m!'b , 






\ a 


-) 








m" 



b — a , , , 
> - A . (56) 



Hence, it follows that there exists no zero point of the function g{x, cj)) in the variable x 
within the interval (6 + 1, 6 + 1 — A). Thus for any a/b G i^Jv-i with 5 < 0, we are able 
to establish that p(6 + 1) = 0. 
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